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Abstract 

With the use of the symmetry-unrestricted cranked Skyrme-Hartree- 
Fock method in the three-dimensional coordinate-mesh representation, 
we have carried out a systematic theoretical search for the superde- 
formed and hyperdeformed rotational bands in the mass A=30-50 re- 
gion. Along the N=Z line, we have found superdeformed solutions in 
32g^ 36 ^j.^ ^°Ca, '^'^Ti, and hyperdeformed solutions in ^^Ar, "^^Ca, 
44rp-^ ^*Cr. The superdeformed band in ^^Ca is found to be extremely 
soft against both the axially symmteric (I30) and asymmetric (1^31) oc- 
tupole deformations. An interesting role of symmetry breaking in the 
mean field is pointed out. 
PACS: 21.60-n; 21.60.Jz; 27.30.+t 
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1 Introduction 



Nowadays, about two hundreds superdeformed (SD) rotational bands are iden- 
tified in various mass (^4=60, 80, 130, 150, 190) regions Every regions 
of superdeformation have their own characteristics so that we can significantly 
enlarge and deepen our understanding of nuclear structure by systematically 
investigating similarities and differences among the SD bands in different mass 
region. For the mass 74=30-50 region, although the doubly magic SD band in 
^^S, which has been expected quite a long time |17|-[15|, has not yet been ob- 



served and remains as a great challenge , quite recently, beautiful rotational 
spectra associated with the SD bands have been observed up to high spin in 
neighboring N=Z nuclei; ^^Ar, ^°Ca, and ^^Ti. In ^^Ar the SD band has 
been identified up to its termination at = 16"*" |T^-|18[. The SD band in 



the spherical magic nucleus ^'^Ca is built on the well known 8p-8h excited 
0"*" states at 5.213 MeV and the rotational spectra have been observed up to 

= 16^ [|19|. In *^^Ti a rotational spectrum associated with the excited 
0+ state at 1.905 MeV has been observed up to I"" = 12+ This rota- 

tional band may also be regarded as belonging to a family of the SD band 
configurations. The fact that rotational bands built on excited 0+ states are 
systematically observed is a quite important, unique feature of the SD bands 
in the ^'^Ca region, as the low angular momentum portions of the SD bands 
in heavier mass regions are unknown in almost all cases. 

In nuclei along the N=Z line, effects of deformed shell structures of pro- 
tons and neutrons act coherently and rich possibilities arise for coexistence 
and competition of different shapes. Thus, we shall be able to learn details of 
deformed shell structure and microscopic mechanism of shape coexistence by 
a systematic study of high-spin yrast structure in the sequence of N=Z nuclei. 
Especially, yrast spectroscopy of nuclei in the ^4=30-50 region, being relatively 
light compared to other regions of SD nuclei, is expected to provide detailed 
information about the roles of individual deformed single-particle orbits re- 
sponsible for the emergence of the SD bands. 

In this paper, as a continuation of the previous work on '^^S |21||, we 
carry out a systematic theoretical search for SD and more elongated hyperde- 
formed (HD) rotational bands in N=Z nuclei from ^^S to ^'^Cr by means 
of the symmetry-unrestricted, cranked Skyrme-Hartree-Fock (SHE) method. 
In Ref. PT], a new computer code was constructed for the cranked SHE cal- 



culation based on the three-dimensional (3D) Cartesian-mesh representation, 
which provides a powerful tool for exploring exotic shapes (breaking both axial 
and reflection symmetries in the intrinsic states) at high spin. The algorithm 
of this code for numerical calculation is basically the same as in Refs. P^-^T 



except that various restrictions on spatial symmetries are completely removed. 
Namely, we do not impose parity and signature symmetries on intrinsic wave 
functions. Hence we call this version of the cranked SHE method "symmetry- 
unrestricted" one. Eor the development of selfconsistent mean-fleld models for 
nuclear structure, we quote Refs. M, p2| and p3l, in which various kinds of 
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mean-field theory, including Hartree-Fock (HF) calculations with finite-range 
Gogny interactions @] and relativistic mean-field approaches are thor- 
oughly reviewd. We also mention that spontaneous symmetry breaking in 
rotating nuclei is reviewd in [0^ . 

In fact, SD and HD solutions of the SHF equations we report in this paper 
preserve the refiection symmetries with respect to the {x,y), {y,z) and {z,x) 
planes, so that the symmetry-unrestricted calculation gives identical results 
with those evaluated by imposing such symmtries. The symmetry-unrestricted 
calculation, however, enables us to examine stabilities of the SD and HD states 
against such refiection-symmetry breaking degrees of freedom like octupole 
deformations. In addition, we shall show that the symmetry breaking play a 
quite interesting role in the crossing region between different configurations 
away from the local minima in the deformation parameter space. 

This paper is arranged as follows: In Section 2, a brief account of the 
cranked SHF method is given. In Section 3, results of calculation for deforma- 
tion energy curves and the SD and HD rotational bands in nuclei from ^^S to 
"^^Cr are systematically presented. Here, special attention will be paied to the 
properties of the SD bands at their high spin limits and the crossover to the 
HD bands with increasing angular momentum. In Section 4, an interesting role 
of symmetry breaking in the mean field will be pointed out in connection with 
configuration rearrangement mechanism. We shall further make a detailed 
analysis of the SD band of ^''Ca and show that it is extremely soft against 
both the axially symmteric (I30) and asymmetric (I31) octupole deformations. 
Main results of this paper are summarized in Section 5. 

A preliminary version of this work was reported in RV, . 



2 Cranked SHF calculation 

The cranked HF equation for a system uniformly rotating about the x-axis is 
given by 

6 < H - uj.otJx >= 0, (1) 

where tUrot and Jx mean the rotational frequency and the x-component of 
angular momentum, and the bracket denotes the expectation value with re- 
spect to a Slater determinantal state. We solve the cranked HF equation for 
a Hamiltonian of the Skyrme type by means of the imaginary-time evolution 
technique |^ in the 3D Cartesian-mesh representation. We adopt the stan- 



dard algorithm p2| , p5H 27|| in the numerical calculation, but completely remove 
various restrictions on spatial symmetries. 

When we allow for the simultaneous breaking of both refiection and axial 
symmetries, it is crucial to accurately fulfill the center-of-mass condition 

AAA 

<Y.^i >=< E Vi >=< E >= 0' (2) 

i=l i=l i=l 
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and the principal- axis condition 



AAA 

< XiHi >=< HiZi >=< Y ZiXi >= 0. (3) 

i=l i=l i=l 

For this purpose we use the constrained HF procedure with quadratic con- 
straints [^]. Thus, we replace the "Routhian" R = H — uJrotJx in Eq- (1) 
with 

3 A 3 A 

R' = R-Yt^k < Y(^k)i >^ - X! f^k,k' < >^ • (4) 

k=l i=l k<k' i=l 

In numerical calculations, we confirmed that the constraints (2) and (3) are 
fulfilled to the order O(10~^^) with values of the parameters /i^ ~ O(IO^) and 
fJ'k,k' ~ 0(1). We solved these equations inside the sphere with radius -R=10 fm 
and mesh size h=l fm, starting with various initial configurations. We note 
that the accuracy for evaluating deformation energies with this mesh size was 
carefully checked by Tajima et al. |]2^,|2^ (see also Ref. [^]) and was found to 



be quite satisfactory. The 9-point formula was used as the difference formula 
for the Laplacian operator. As usual, the angular momentum is evaluated as 
Ih =< Jx >■ For the Skyrme interaction, we adopt the widely used three 



versions; Sill |T1, SkM* H and SLy4 |3 . 

In addition to the symmetry-unrestricted cranked SHF calculation ex- 
plained above, we also carry out, for comparison sake, symmetry-restricted 
calculations imposing refiection symmetries about the (x, y)-, {y, z)- and {z, x)- 
planes. The computational algorithm for this restricted version of the cranked 



SHF calculation is basically the same as in , but we have constructed a new 
computer code for this purpose. Below we call these symmetry-unrestricted 
and -restricted cranked SHF versions "unrestricted" and "restricted" ones, re- 
spectively. Comparison between results obtained by unrestricted and restricted 
calculations carried out independently serves as a check of numerical results to 
be presented below. Physical significance of this comparison is, however, that 
we can, in this way, clearly identify effects of symmetry-breaking in the mean 
field. We shall indeed find an interesting symmetry breaking effect in the next 
section. 

Solutions of the cranked SHF equation give minima in the deformation 
energy surface. In order to explore the deformation energy surface around 
these minima and draw deformation energy curves as functions of deformation 
parameters, we carry out the constrained HF procedure with quadratic con- 
straints ||39|. Namely, in addition to the constraints to fulfil the center-of mass 



and principal-axis conditions mentioned above, we also introduce constraints 
involving relevant mass-multipole moment operators and solve resulting con- 
strained HF equations. 

As measures of the deformation, we calculate the mass-multipole moments, 

47r f 

aim = J r^Xim (n) p (r) dr, (m = -Z, ■ ■ ■ , /) (5) 
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where p{r) is the density, R — < Z^^i r*f > /S^l, and Xi^ are real bases of 
the spherical harmonics, 



Xio = 1/0, (6) 
1 



Xi\m\ = ^(F;_|r„| + (7) 



^i-H = ;^(^lmi - (8) 

Here the quantization axis is chosen as the largest (smallest) principal axis for 

prolate (oblate) solutions. We then define the quadrupole deformation param- 
eter /32, the triaxial deformation parameter 7, and the octupole deformation 
parameters [3^ and [3zrn by 



Oi2o = 1^2 COS 7, = 1^2 sin 7, (9) 



1/2 

, /?3m= («L + «3-m)''" (m = 0, 1, 2, 3) . (10) 




For convenience, we also use the familiar notation — for oblate shapes 
with (/32,7 = 60°). 

3 Results of calculation 

3.1 Deformation energy curves 

Figures 1-3 show deformation energy curves evaluated at / = by means 
of the constrained HF procedure with the quadratic constraint on the mass- 
quadrupole moment. The SIII, SkM*, and SLy4 versions of the Skyrme in- 
teraction are used in Figs. 1, 2, and 3, respectively. Solid lines with and 
without filled circles in these figures represent results of unrestricted and re- 
stricted calculations, respectively. Let us focus our attention to the region of 
large quadrupole deformation ^2- In both cases, we otain local minima cor- 
responding to the SD states for ^^S, ^^Ar, ^°Ca and ^^Ti in the region 
0.4 < (52 < 0.8. (The local minimum in ^^Ti is triaxial, as shown in Fig. 10 
below, i.e., it is situated away from the 7 = section of the deformation en- 
ergy surface, so that it is not clearly seen in Figs. 1-3.) The local minima in 
^^S and ^^Ar involve four particles (two protons and two neutrons) in the 
fp shell, while those in ^°Ca and ^^Ti involve eight particles (four protons 
and four neutrons). These local minima respectively correspond to the 4p-12h, 
4p-8h, 8p-8h and 8p-4h configurations with respect to the doubly closed shell 
of ^°Ca, and their properties have been discussed from various point of view; 
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see Refs. [0-0 for ^^S, Refs. |T|-|T8|,g for ^^Ar, Refs. [|1|,|5|-|3] for ^^Ca, 
and Refs. [10,11,11] for ^^Ti. 

In addition to these SD minima, we also obtain local minima in the region 
(^2 > 0.8 for ^°Ca, ^^Ti and ^^Cr. These minima involve additional four parti- 
cles (two protons and two neutrons) in the single-particle levels that reduces to 
the gg/2 levels in the spherical limit. Somewhat loosely we call these local min- 
ima "hyperdeformed." The HD solution in ^°Ca corresponds to the 12p-12h 
configuration. For ^^Ti, we obtain two HD solutions which correspond to the 
12p-8h and 16p-12h configurations. The HD solution in ^^Cr corresponds to 
the 16p-8h configuration. These HD solutions well agree with those previously 
obtained in the SHF calculation by Zheng, Zamick and Berdichevsky [^]. We 
also mention that the 12p-12h configuration in ^°Ca and the 16p-12h con- 
figuration in '^^Ti agree with those obtained by the macroscopic-microscopic 
model calculation by Leander and Larsson 

In Figs. 1-3 and in the following, the above SD and HD configurations are 
denoted by (or (/p)"*?™), where n and m indicate the numbers of nucleons 
occupying the fi/2 shell (or the fp shell) and the gQ/2 shell, respectively. 

As seen in Figs. 1-3, these SD and HD minima are obtained for all cal- 
culations with the use of the SHI, SkM*, and SLy4 interactions. These local 
minima preserve the reflection symmetries so that the results of restricted and 
unrestricted calculations are the same. On the other hand, we also find a case 
where the two calculations give different results: We obtain a HD minimum 
with /?2 — 0.8 for ^^Ar in the restricted calculations. This minimum involves 
eight particles (four protons and four neutrons) in the fp shell and correspond 
to the 8p-12h configuration, but it disappears in the unrestricted calculations 
and its remnant remains as a shoulder of the deformation energy curve. 

Although the restricted and unrestricted calculations give identical results 
for the SD and HD local minima except for the HD solution for ^^Ar, they 
show different behaviors in regions away from the local minima: In Figs. 1-3, 
we see that the deformation energy curves obtained by the unrestricted calcu- 
lations always join different local minima smoothly. On the other hand, in the 
restricted calculations, segments of the deformation energy curves associated 
with different local minima sharply cross each other in some situations, while 
they are smoothly joined in other situations. Closely examining the configura- 
tions involved, we notice that the sharp crossings occur between configurations 
having different numbers of particles excited into the fp shell. This point will 
be further elaborated in the subsequent section. 

In Figs. 1-3 there are a number of local minima in the region of smaller 
values of (32- We shall not discuss on these local minima in this paper, since the 
pairing correlations not taken into account here are expected to be important 
for these. 

The single-particle level schemes at the SD and HD local minima men- 
tioned above are displayed in Figs. 4-6 at positions of their quadrupole defor- 
mation parameters f32- One may notice that these equilibrium deformations 
correspond to values of f32 slightly smaller than those where the amount of 
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energy spacings between the highest occupied and lowest unoccupied single- 
particles levels become local maxima. This is a characteristic common to 
three versions of the Skyrme interactions used, and in accord with the expec- 
tation that equilibrium deformations are determined by the sum of microscopic 
(shell-structure) and macroscopic (liquid-drop) energies; the latter shifts the 
equilibrium values of (32 to slightly smaller ones. 



3.2 SD and HD rotational bands 

Let us focus our attention to the SD and HD local minima shown in Figs. 1- 
3, and investigate properties of the rotational bands built on them. Figures 
7-9 show excitation energies, as functions of angular momentum, of the SD 
and HD rotational bands calculated with the use of the SHI, SkM*, and SLy4 
interactions, respectively. These rotational bands are obtained by cranking 
each SHF solution (the SD and HD local minima in Figs. 1-3) and following 
the same configuration with increasing value of Wrot until the point where we 
cannot clearly identify the continuation of the same configuration any more. 
We note that, in "^^Ti, two HD bands associated with the f^g"^ and {fpY'^g^ 
configurations cross at / = 30 — 34, and the latter becomes the yrast for higher 
spin. (This band continues beyond J = 40 where the figure is cut.) In addition 
to the SD and HD bands built on the = 0+ band-head states, we have found 
a HD band in ^^Ar, which does not exist at / = and emerges at / ~ 16 
due to the rotation alignment of the gg/2 orbit. This HD band is denoted by 
f^g'^ and is included in Figs. 7-9. A similar configuration, f^g'^, was found for 



^^S in our previous calculation I^T] and called "HD-like." This and analogous 
configurations in nuclei other than '^^Ar are not illustrated in Figs. 7-9 in 
order not to make the figure too complicated (drawing complete yrast spectra 
of individual nuclei is not the major purpose of these figures). 

As is well known, according to the deformed harmonic-oscillator potential 
model, = Z = 18 and 24 are magic numbers associated with the HD shell 
structure with axis ratio 3:1. These HD states respectively correspond to 
the f'^g'^ and f^'^g'^h'^ configurations in our notation, where h denotes the level 
associated with the /iii/2 shell. Microscopic structures of the HD solutions 
under discussion are apparently different from these, however. We also mention 
that the possible existence of HD rotational bands at high spin in ^^Ar and 
^^Cr have been discussed in Refs. [^,|5^ from the viewpoint of the cranked 



cluster model. The relationship between our solutions and their solutions 
associated with cluster structure is not clear. 

Calculated quadrupole deformation parameters (/32,7) of all bands men- 
tioned above and their variations are displayed in Fig. 10. The rotational 
frequency dependence of the single-particle energy levels (Routhian) is illus- 
trated in Fig. 11, taking the SD band in ^°Ca as a representative case. The 
excitation energies of the SD and HD bands obtained by using different ver- 
sions (SHI, SkM*, SLy4) of the Skyrme interaction are compared in Fig. 12 
with the experimental data []TB|,^,^. 
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Examining these figures, we see that, aside from quantitative details and 
some subtle points to be discussed below, the results obtained by using different 
versions of the Skyrme interaction are similar. This implies that the basic 
properties of the SD and HD bands under discussion are not sensitive to the 
details of the effective interaction. 

As mentioned in the introduction, one of the unique features of the SD 
bands in the ^°Ca region is the possibility to observe the SD rotational level 
structure from the I'" = O"*" band heads up to the maximum angular momenta 
allowed for the many-particle-many-hole configurations characterizing the in- 
ternal structures of these bands. In fact, such a "SD band termination" has 
been observed at / = 16 in ^^Ar and well described by calculations in terms 



of the j-j coupling shell model, the cranked Nilsson-Strutinsky model p!6|-p!8 



and the projected shell model [44|. On the other hand, for ^°Ca and ^^Ti, 
it is not clear whether or not the SD band continues beyond the highest spin 
states observed up to now (the 16"*" state in "^^Ca and the 12"'" state in 



^^Ti [^]) and, quite recently, their properties, from the 0"*" band- heads to such 
high-spin regions, have been discussed in terms of the spherical shell model in 



Ref. [51 1 for Ca and in Ref. [20 for Ti. In our calculation, except for 



the case of using the SLy4 interaction, the band termination phenomenon in 
^^Ar is reproduced; the shape becomes triaxial and evolves toward the oblate 
shape, although the oblate limit is not reached. In the cases of ^°Ca, the 
shape is slightly triaxial with 7 = 6°-9°(8°-9°) and the SD band terminates at 
/ ~ 24 in the calculation with the use of the SIII (SkM*) interaction. In the 
case of "^^Ti, the shape is more triaxial with 7 = 18°-25° and 13°-19°, and 
the SD bands terminates at J ~ 12 and 16 for the SIII and SkM* interactions, 
respectively. Thus, the band termination properties appear quite sensitive to 
the details of the effective interaction. Concerning the SD band termination 
in ^°Ca and ^^Ti, the results obtained with the use of the SIII and SkM* 
interactions would be more reliable than that with SLy4, in view of the above 
discussion for ^^Ar. In any case, it would be very interesting to explore higher 
spin members of the SD rotational bands in ^°Ca and ^^Ti in order to 
understand the terminating properties of the SD bands at high spin limits. 

As is clear from the comparison with experimental data in Fig. 12, the mo- 
ments of inertia for the SD band are somewhat overestimated in the present 
calculation. To investigate a possible cause of this, we plan to take into ac- 
count the pairing correlations by means of the cranked Skyrme-Hartree-Fock- 
Bogoliubov code constructed in Ref. [Q. One also notice that the excitation 
energy of the SD band-head state in ^°Ca is overestimated. It will decrease 
if the zero-point rotational energy correction, —jj < >, is taken into ac- 
count (see Ref. [4£] for numerical examples). Although the calculation of this 
correction is rather easy, we need to evaluate, for consistency, also the zero- 
point vibrational energy corrections [0, and this is not an easy task. We 
therefore defer this task for a future publication. Inclusion of these correla- 
tions is expected to improve agreement with the experimental data. 



8 



4 Discussions 



4.1 A role of symmetry breaking 

Let us now discuss on the significance of tlie reflection symmetry breaking in 
the mean field. As noticed in Figs. 1-3, the crossings between configurations in- 
volving different numbers of particles in the fp shell are sharp in the restricted 
calculation, while we always obtain smooth configuration rearrangements in 
the unrestricted calculations. The reason for this different behavior between 
the unrestricted and restricted calculations is rather easy to understand: When 
the parity symmetry is imposed, there is no way, within the mean-field approx- 
imation, to mix configurations having different number of particles in the fp 
shell. In contrast, smooth crossover between these different configurations is 
possible via mixing between positive- and negative-parity single-particle lev- 
els, when such a symmetry restriction is removed. Let us examine this idea in 
more detail. In Figs. 13-15 octupole deformation parameters /?3 of the lowest 
energy states for given values of (32 are shown. They are obtained by the unre- 
stricted SHF calculations and plotted as functions of P2 in the lower portion of 
each panel. We see that (3^ are zero near the local minima in the deformation 
energy surface, but rise in the crossing region between configurations involving 
different number of particles in the fp shell. This means that the configuration 
rearrangements in fact take place through paths in the deformation space that 
break the reflection symmetry. The importance of allowing the mean field for 
breaking symmetries in the process of configuration rearrangements was pre- 
viously emphasized by Negele in their calculations for spontaneous fission 
of ^^S by means of the imaginary time tunneling method. 

In connection with the above finding, it may be appropriate to point out 
another situation in which a symmetry breaking in the mean field plays an 
important role. For both the restricted and unrestricted calculations, we have 
obtained smooth crossover between the SD and HD configurations in ^°Ca and 
^^Ti (see Figs. 1-3). Since four particles are further excited into the g shell 
in the HD configurations, the smooth configuration rearrangement becomes 
possible by means of the mixing between the down-sloping levels steming 
from the gg/2 shell (its asymptotic quantum number is [440] |) and the up- 
sloping levels stemming from the sd shell ([202]| and [200] | in the cases of 
^°Ca and ^^Ti, respectively). The mixing between these single-particle levels 
takes place through the hexadecapole components of the mean field, and we 
need to break the axial symmetry to mix them in the case of ^°Ca. The cal- 
culations called "restricted" in this paper allow the axial symmetry breaking, 
so that the smooth rearrangement from the SD to HD configurations is pos- 
sible also in ^°Ca. A very careful computation is required, however, in order 
to detect these mixing effects, since the interaction between the down-sloping 
and up-sloping levels is extremely weak. 

In Figs. 13-15, one may notice that (3^ take non-zero values also in some sit- 
uations other than the crossing regions. Such situations occur in some regions 
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of the deformation energy surface where it becomes very soft with respect to 
the reflection-asymmetric degrees of freedom. In the next subsection, we in- 
vestigate this point in detail taking the SD solution in ^°Ca as an especially 
interesting example. 

4.2 Octupole softness of the SD hand in ^^Ca 

Let us examine stabilities of the SD local minimum in ^°Ca against octupole 
deformations. Figure 16 shows deformation energy curves as functions of the 
octupole deformation parameters P^mijn = 0, 1,2,3) for fixed quadrupole de- 
formation parameters at and near the SD minimum of ^''Ca, calculated by 
means of the constrained HF procedure with the use of the SIII, SkM*, and 
SLy4 interactions. We immediately notice that the SD state is extremely soft 
with respect to the P^o and Psi deformations, irrespective of the Skyrme in- 
teractions used. Although it is barely stable with respect to these directions 
(see curves for P2 = 0.6), an instability toward the P^i deformation occurs as 
soon as one goes away from the local minimum point (see curves for (32 = 0.5). 
In fact, the deformation energy surface is found to be almost flat for a com- 
bination of the jJ^o and ^31 deformations already at the SD local minimum. 
Thus we need to take into account the octupole shape fluctuations for a better 
description of the SD rotational band in ^°Ca. It will be a very interesting 
subject to search for negative-parity rotational bands associated with octupole 
shape fluctuation modes built on the SD yrast band. We plan to make such 
a study in future. Quite recently, the octupole instability of the SD band in 
^°Ca has been suggested also by Kanada-En'yo [ pO[1 . 

5 Conclusions 

With the use of the symmetry-unrestricted cranked SHF method in the 3D 
coordinate-mesh representation, we have carried out a systematic theoretical 
search for the SD and HD rotational bands in the N=Z nuclei from ^^S to 
^^Cr. We have found the SD solutions in ^^S, ^^Ai, ^°Ca, ^^Ti, the HD 
solutions in ^^Ar, ^°Ca, ^^Ti, ^'^Cr, and we have carried out a systematic 
analysis of their properties at high spin. 

It is explicitly shown that the crossover between configurations involving 
different number of particles in the fp shell takes place via a reflection-symmery 
breaking path in the deformation space. 

Particular attention has been paid to the recently discovered SD band in 
^°Ca, and we have found that the SD band in ^'^Ca is extremely soft against 
both the axially symmetric (I30) and asymmetric (Vsi) octupole deformations. 
Thus, it will be very interesting to search for negative-parity rotational bands 
associated with octupole shape vibrational excitations built on the SD yrast 
band. 
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Figure 1: Deformation energy curves as functions of the quadrupole deformation 

(32 calculated at I = by means of the constrained SHF procedure with the SIII 
interaction. The axial-asymmetry parameter 7 is constrained to be zero. The curves 
for different nuclei are shifted by 20 MeV to accommodate them in a single plot. Solid 
lines with and without filled circles represent the results obtained by the unrestricted 
and restricted versions, respectively (see the text). The notations f^g"^ and {fp)^g"^ 
indicate the configurations in which the fj/2 shell (fp shell) and the gg/2 shell are 
respectively occupied by n and m nucleons. 
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Figure 2: The same as Fig.l but for the SkM* interaction. 
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Figure 3: The same as Fig.l but for the SLy4 interaction. 
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Figure 4: Single-particle energy spectra for neutrons at various SD and HD local 
minima, calculated at / = with the SIII interaction. For every local minima, 
levels arc drawn at positions corresponding to their equilibrium values of (32- The 
levels are occupied up to the Fermi levels; on top of them, neutron numbers of the 
system are indicated. The levels associated with the fi/2 and gg/2 shells are drawn 
by broken and thick solid lines, respectively. Single-particle spectra for protons are 
almost the same as those for neutrons. 
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Figure 5: The same as Fig.4 but for the SkM* interaction. 
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Figure 6: The same as Fig.4 but for the SLy4 interaction. 











Ar-.fV 




Ti -(.fp)'V 






















Ti-.fV 




Ca-f 




Ar-f 







5 10 15 20 25 30 35 40 
Angular Momentum l[h] 



Figure 7: Excitation energy vs. angular-momentum plot for the SD and HD ro- 
tational bands obtained by the cranked SHF calculations with the use of the SIII 
interaction. Their configurations are indicated by the same notations as in Fig. 1. 
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Figure 8: The same as Fig.7 but for the SkM* interaction. 
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Figure 9: The same as Fig.7 but for the SLy4 interaction. 
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Figure 10: Quadrupolc deformation parameters {P2i7) and tlieir variations for the 
SD and HD rotational bands in ^'^S, ^^Ar, ^°Ca, ^"^Ti and "^^Cr. The top, middle 
and bottom panels show the shape evolution in the {P2, 7) plane, evaluated with the 
use of the SIII, SkM*, and SLy4 interactions, respectively. Notations and ranges 
of / for individual bands are the same as those shown in Figs. 7-9. Values of P2 
decrease with increasing I, except for the (/p)^^^^ configuration in ^"^Ti (where P2 
increases with increasing I). 
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Figure 11: Neutron single-particle energy diagram (Routhian) for the SD band in 
^°Ca, calculated with the use of the SIII interaction and plotted as functions of 
rotational frequency Wrot- The levels associated with the /7/2 and gg/2 shells are 
drawn by thick-broken and long-dashed lines, respectively. 
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Figure 12: Comparison between the excitation energies of the SD and HD bands 
in ^^Ar, ^'^Ca, ^^Ti and ^^Cr, calculated by using different versions of the Skyrme 
interaction and the experimental data ( [16| for ^^Ar, for ^'^Ca, and |^0| for 
^^Ti). The data are shown by filled circles and the results with the SIII, SkM* and 
SLy4 interactions are drawn by solid, dashed and dashed-dotted lines, respectively. 




Figure 13: Octupole deformation paramaters /^s of the lowest energy states for 
given values of /32, obtained by the unrestricted SHF calculation at / = with the 
use of the SIII interaction. Their values are plotted as functions of (32 in the lower 
portion of each panel. To show that their values increase at crossing regions between 
configurations involving different number of particles in the fp shell, the deformation 
energy curves are also displayed in the upper portion of each panel. The latter are 
the same as those presented in Fig. 1. 
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Figure 14: The same as Fig. 13 but for the SkM* interaction. 




Figure 15: The same as Fig. 13 but for the SLy4 interaction. 
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Figure 16: Left-hand side: Deformation energy curves (measured from energies at 
P3 = 0) as functions of the octupole deformation parameters Psmini = 0,1,2,3), 
calculated for ^'^Ca by means of the constrained HF procedure. The quadrupole 
deformation parameters are fixed at /32 = 0.5 and 7 = 0. One of the Psmii^ = 
0,1,2,3) is varied while the other /^sm's are fixed to zero. Right-hand side: The 
same as the left-hand side, except that the quadrupole deformation parameters are 
fixed at P2 = 0.6 and 7 = 0. Results of calculation with the use of the SIII, 
SkM*, and SLy4 interactions are displayed in the upper, middle and lower panels, 
respectively. 



21 



